Quantum theory of light and noise polarization in nonlinear optics 
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We present a consistent quantum theory of the electromagnetic field in nonlinearly responding 
causal media, with special emphasis on x^^^ media. Starting from QED in linearly responding 
causal media, we develop a method to construct the nonlinear Hamiltonian expressed in terms of 
the complex nonlinear susceptibility in a quantum mechanically consistent way. In particular we 
show that the method yields the nonlinear noise polarization, which together with the linear one is 
responsible for intrinsic quantum decoherence. 
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Recent advances in quantum information technologies 
have been the main driving forces behind the desire to 
build parametric down- conversion sources of entangled 
photon pairs (in the low-intensity limit) or 
two-mode squeezed states (in the high-intensity limit) 
with high fidelity. It is known that single-photon states of 
non-unit efficiency as produced by heralded single-photon 
sources using parametric down-conversion, cannot be pu- 
rified by using linear optical elements andphoto detec- 
tion to yield states with higher efficiency |3, 1^ • That in 
turn means that post-processing of single-photon sources 
is impossible and the sources themselves have to be im- 
proved. In order to achieve the maximally possible purity 
of heralded single-photon states or correlated (entangled) 
twin-beam photons it is therefore necessary to investigate 
the theoretical limits nature imposes on us. 

An important step in this direction is to provide a 
quantum theory of light that takes into account nonlin- 
ear processes such as parametric down-conversion, and 
at the same time decoherence mechanisms due to un- 
avoidable absorption losses of the nonlinear material the 
light interacts with. The theory of quantized electro- 
magnetic fields in linearly and causally responding ma- 
terials (with the linear response function satisfying the 
Kramers-Kronig relations) is well established (see, e.g., 
Refs. p. [lol[lll |). It has been known for some time that 
analogous Kramers-Kronig relations do also hold for non- 
linear susceptibilities Hence, it will be interesting 
to see how these causal relations appear in a nonlinear 
quantum theory. 

Previous work on electromagnetic field quantization in 
nonlinear materials have focused on strictly lossless ma- 
terials where Lagrangian methods and mode decompo- 
sitions apply [U Im lisl llH . A first attempt to in- 
clude in the field quantization both linear and nonlin- 
ear losses was made in Ref. for Kerr media, by ex- 
tending the linear harmonic-oscillator model used in the 
Huttner-Barnett quantization scheme to a nonlinear 



one. A consistent approach that includes — for given non- 
linear susceptibility — absorption and dispersion has not 
yet been formulated within the frame of (macroscopic) 
QED. 

In this article we will exemplify, on focusing on x^^-* 
media, how to consistently quantize the electromagnetic 
field in the presence of nonlinearly responding causal ma- 
terials. This theory provides the starting point for further 
investigations of theoretical limits to the performance of 
nonlinear optical elements as sources of nonclassical light. 
Starting from the nonlinear Hamiltonian expressed in 
terms of the canonically conjugated variables as used in 
QED in linear causal media, we first express the nonlinear 
polarization field in terms of these variables as well. This 
is compared with the classical nonlinear response which 
enables us to identify the nonlinear noise contributions. 

We begin with recapitulating the quantization scheme 
for the electromagnetic field in the presence of a linearly 
(and locally) responding causal dielectric medium of per- 
mittivity e(r, w) =£'(r, w) -I- ie"(r, w) lOj. In this case the 
Hamiltonian is bilinear. 
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with the annihilation and creation operators (r, oj) and 
JI{y,lo), respectively, playing the role of the canonically 
conjugated dynamical variables which are attributed to 
collective excitations of the electromagnetic field and 
the dispersing and absorbing dielectric matter and obey- 
ing the bosonic commutation rules [/^(r, w), /^(r',^')] = 
6ijS{LU — Lu')S{r — r'). By expressing the electromagnetic 
field in terms of the dynamical variables, the electric field, 
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for example, reads 



E(r) = / dcjE(r,w) +h.c., (2) 



E(r,c^) - 




(3) 

it can then be shown that Maxwell's equations, in partic- 
ular Faraday's and Ampere's laws, hold. In Eq. the 
(dyadic) Green function G(r, s, w) is the unique funda- 
mental solution of the inhomogeneous Helmholtz equa- 
tion 

V X V X G(r,s, w)-— e(r, w)G(r,s, = S{r-s)I (4) 

and contains all relevant information about the material 
properties and the geometry of the system. 

Equation ^ together with Eq. Q may be regarded as 
a generalization of the ordinary mode expansion, with the 
role of the mode operators being taken on, in a sense, by 
the fi{r,uj,t) and fj{r,oj,t). In summary, (i) the Hamil- 
tonian generates, using the representation according 
to Eqs. (0) and Q, the correct (macroscopic) Maxwell 
equations, (ii) the fundamental QED equal-time commu- 
tation relations are preserved, and (iii) the fluctuation- 
dissipation theorem takes its standard form. Moreover, 
the Hamiltonian represents an energy stored in the 
system composed of the electromagnetic field and absorb- 
ing matter. 

To turn over to the nonlinear media, let us first fix 
some notation. From now on we will abbreviate spatial 
and frequency variables (rk,ujk) by their label fe, e.g., 
l = (ri,LJi) and write J dk = J d^rt J dtuk- In the latter 
integrals, the spatial integration extends over all space. 
The frequency integral, which we initially will assume to 
range over all positive frequencies, will be restricted later 
on. On recalling the physical meaning of the dynamical 
variables fi{r,uj) and fj{r,uj), the most general normal- 
order form of the nonlinear interaction energy that cor- 
responds to a x*-^-* medium reads 

Hnl = y"dld2d3a,(,fc)(l,2,3)//(l)/,(2)/fc(3)+h.c.. 

(5) 

The unknown tensor function otujk) (1, 2, 3), which has to 
be symmetrized over its last two indices to avoid double- 
counting, has to be determined from constraints imposed 
by generally accepted relations. 

We first note that Faraday's law can be written as 

V X E(r) = -B(r) = -1 [B(r), Hl + H^l] ■ (6) 
in 

Both the electric and the magnetic induction fields are 
pure electromagnetic fields without being related to the 
material degrees of freedom and hence their equal-time 



commutation relations are as in vacuum QED. To be 
more specific, we may assume that the functional form 
of these fields in terms of the dynamical variables f(r, oj) 
and fj{r,uj) is (in close analogy to the case of ordinary 
mode expansion) the same as in the linear theory. From 
here it immediately follows that 

[B(r),HjvL] =0. (7) 

Using Faraday's law, we rewrite Ampere's law as 

V X V X E(r) = -MoDL(r) - MoPjVL(r), (8) 

where we have split up the dielectric displacement field 
D(r) into the linear part Di(r) and the nonlinear po- 
larization PjvL(r)- Employing Heisenberg's equation of 
motion, we may rewrite Eq. JSJ as 

Vx VxE(r)-0[[DL(r),i7i],i/i] 

^ ^{[[nLir),HL],HNL] + [[T>Lir),HNL],HL] 
+ [[PNLir),HL],HL]}, (9) 

where we have kept, for consistency reasons, only terms 
that are in at most first order in the nonlinear coupling 
coefficient ai(jfe) (1, 2, 3). The Ihs of Eq. is zero by 
the definition of the linear displacement field. Note that 
the time dependence is carried by the time-dependent 
dynamical variables fi{r,u!,t) and f}{Y,LO,t). The first 
term on the rhs of Eq. Q vanishes by virtue of the 
constraint 0. To see this, one has to express the lin- 
ear displacement and the magnetic induction fields in 
terms of the dynamical variables, leading to [DL(r), Hl] 

— {ih/ X B(r), and application of Eq. Q leads to 
the quoted result. Hence, we are left with a relation 
between double commutators of the linear displacement 
and nonlinear polarization fields with the linear and non- 
linear parts of the Hamiltonian, [[T)L{r), Hnl], Hl] = 

— [[Pnl{j^), Hl], Hl]- A particular solution is certainly 

[T>L{r),HNL]-^ -[PNLir),HL]. (10) 

The general solution would additionally include commu- 
tants with the linear Hamiltonian Hl- These terms 
must be functionals of the number (density) operator 
f^(r,a;) • i{r,uj)- However, we can assume that linear 
functionals of this type are already included in the partic- 
ular solution (|10|l as they lead to bilinear forms in the dy- 
namical variables. On the other hand, quartic and higher 
functionals have to be excluded to ensure that Pjvi(r) 
stays bilinear which guarantees consistency within the 
approximations made. 

The expression on the rhs of Eq. 1)10(1 is nothing but 
the Liouvillian Cl generated by the linear Hamiltonian 
Hl acting on the nonlinear polarization field. Therefore, 
Eq. l(Tm) can be solved for PjvL(r) to yield 

PNL{r) = ^^Cl'[nL{r),HNL]- (11) 
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At this point we recall that according to 

= eoe(r,c.)E(r,c.)+Pf^(r,c.) (12) 



the linear displacement field 



DlW = / (iu;Di(r,w) +h.c. 
Jq 



(13) 



consists of a reactive part related to the electric field and 
a noise part 



(r) 



dcj P}j^\r, uj) +h.c. 



(14) 



Inserting Eq. lfO|l together with Eq. (fT^ into Eq. (fTT|l . 
we see that the nonlinear polarization also decomposes 
into a reactive part, which can be related to the non- 
linear response, and a noise part, which determines the 
nonlinear noise polarization 



NL 



(r) 



:C-, 



Because of the relation P^^^(r,a;) = iy^heo/Try^e"{r, uj) 

f(r,a'), P^^(r) vanishes if the imaginary part of the lin- 
ear permittivity, s"{rAj), and hence the noise associated 
with it tends to zero |20j. 

The inverse Liouvillian can be calculated using stan- 
dard techniques, and we obtain from Eq. Hll|) 

PjVL(r) = 

lim /dre-'^e-^^^^ [Di,(r),i/jVL]e-^^^^, (16) 



NL 



(15) 



where the real positive number s ensures convergence of 
the integral. In the next step we compute the commuta- 
tor [DL{r), Hnl] and evaluate the integral in Eq. (|16|l . 
First, we evaluate the commutator between the dynami- 
cal variables and the nonlinear Hamiltonian H^l, lead- 
ing to [here, = (s, lu)] 

[/„(0),Hwi] = J d2rf3a,„(,fe)(0,2,3)/,(2)A(3) 

+ 1 (ild2a*(,„)(l,2,0)/](2)/,(l). (17) 

In what follows, we will concentrate on the contribu- 
tion to the nonlinear displacement and polarization that 
comes from terms containing two annihilation opera- 
tors such as /j(2)/fe(3). We will label these contri- 
butions with the superscript ^^~^'> in analogy with the 
standard notation for positive-frequency parts. The in- 
verse Liouvillian of the bilinear combination of annihi- 
lation operators is readily found to be fj{2)fk{3) = 



i/(w2+W3)/j (2)/fe(3). Combined with Eq. ((TB|l we finallv 
obtain for the nonlinear polarization field 



P, 



UJ2 + 



O(0,2,3) 



X ^£(r,c.)Gi™(r,0)/,(2)/,(3)+P^^';+)(r), (18) 
where the noise polarization reads 

X ^(r-s)/,(2)/fc(3). (19) 

In order to make contact with standard notation, let us 
recall the definition of the nonlinear polarization within 
the framework of response theory: 

t 

Pnl,i{^, t) = £o J dTidT2 Xhi„(r, t - n, i - T2) 

X ii;,„(r, Ti)i?„(r, T2) -f P\^li[v, t). (20) 

The first term on the rhs of Eq. (|20l) is the causal response 
well known from nonlinear optics [l^, with x|^„(r, ti, i2) 
being the response function of the x^^-* medium. The 
term P^'j(r, t) is a (yet unknown) nonlinear noise polar- 
ization commonly disregarded in classical nonlinear op- 
tics. In most cases of interest it is sufficient to evaluate 
Eq. 1)20(1 in the slowly-varying amplitude approximation 
in the sense that 



E(r,t) = ^E,(r,f7^,t)e 



~iQ.„t 



(21) 



with the time scale on which the amplitude function 
Ei{v,^l^,t) noticeably changes being long compared 
with VL~^ and the characteristic time of variation of 
X;V„„(r, ^1, ^2) with respect to both ti and t2 (see, e.g., 
the treatment in Ref. 0|). Hence the slowly varying 
field amplitudes can be taken out of the integral at the 
upper integration limit i, and we are left with the Fourier 

transform of x|™„(r, ti, ^2), x|™„(r, t^i, ^^2), which slowly 
varies with uji and UJ2- In this way we derive 



P^++^(r,f^23)-eox[l(r,f^2,f^3) 



X E^{v,^2)En{v,^^) + p\^^i{YM 



NL,l\ 



'■23) 



(22) 



[r223 = f22 + ^^3] , where the time argument t of the ~ quan- 
tities has been omitted for notational convenience. 

The validity of the approximation leading from 
Eq. to Eq. may be regarded as being a prereq- 
uisite for substantiating the effective interaction Hamil- 
tonian (jSJ. At the same time, it suggests further specifi- 
cation of the Hamiltonian as therein the introduction of 
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slowly varying variables is desirable. In view of Eqs. ^ 
and we define, on assuming the Green tensor and 
the linear susceptibility are slowly varying with uj, the 
slowly varying bosonic variables i{r,n^) = (Afij,)"^/^ 
Iaq dLui{v,uj,t)e^^^^ (Aili,, relevant frequency interval 
around fii/), and Eq. © reduces to 

HnL = J d^Sid^S2(f S3 ai(^jk)isi, ^23:^2:^2, S3, fls) 

^/MhMhMh (si , 1]23 ) /, (S2 , 2 ) /fc (S3 , 3 ) + h.c. . 

(23) 

Introducing in Eqs. (|18|l and (|22|l the slowly vary- 
ing variables {(r,n^), from a comparison of the reac- 
tive parts of the nonlinear polarization as given by the 
two equations we derive the following integral equa- 
tion for determining the nonlinear coupling coefficient 
<^i{jk) (si, f^23, S2, ^2', S3, ^^3) in terms of the nonlinear sus- 
ceptibility x[m„(r, ^^2, ^^3): 

J d^s ^e"{s,n23) am(jk) (s, ^^23, S2, fi2; S3, ^3) 

X G;„i(r,s,1723) - ' ^ ^ 



ITTC^ V feo ^^23 



X v/£"(s2,fi2)£"(s3,1^3)xll(r,f^2,r!3) 

X G™j(r,S2,r22)G„fe(r,S3,f^3). (24) 

This equation is of Fredholm type and can be solved by 
inverting the integral kernel on the Ihs of Eq. (I24II . Note 
that the inverse of the Green tensor is just the Hclmholtz 
operator Hij{r,uj) = d^dj — 5ij/V — (w^/c^)e(r, ct;)^^ : 
i?jj(r, a;)Gjfc(r, s, w) = (5ifcJ(r — s). Hence, from Eq. 
it follows that 

2 nl^l 



I — ~ — 

a,(,,)(r,l]23,S2,f^2;S3,l)3) = — 



/e"(s2,r!2)e"(s3,f^3) 



e"(r,f]23) e(r,f^23) 



^23 
Hu{y, ^23) 



X [xli(r,r!2,f^3)G„„(r,S2,r!2)G„fc(r,S3,l}3)]. (25) 

Re-inserting Eq. H25|) into Eq. p9|) eventually yields, 
on recalling Eqs. (0), ®, H21|l and the definition of the 
slowly varying variables f (r, £7^), the following expression 
for the nonlinear noise polarization: 



^NL.l (.r,ii23j 



nl^e{v,n23) 



X HH(r,r!23)[xl;^„(r,r!2,f^3)i?™(r,r!2)i?„(r,173)]. (26) 

To our knowledge, this is the first time a nonlinear 
noise polarization has been derived in the frame of quan- 
tum nonlinear optics. Note that the Helmholtz opera- 
tor acting on the electric field returns the linear noise 
polarization, Hij{Y,uj)E_j{Y,uj)=uj'^ / {eQC^)p}i^} (v,uj) [cf. 



Eq. (|12(l ]. Among other terms, Eq. H26|) contains prod- 
ucts of the electric field and the linear noise polarization. 

In summary, we have presented a consistent quan- 
tum theory of the electromagnetic field in the pres- 
ence of quadratically responding dielectric materials. It 
takes care of the causal nature of the dielectric response 
which implies the existence of a nonlinear noise polar- 
ization. The nonlinear (effective) interaction Hamilto- 
nian H23|) [or equivalently, Eq. jSJ in the slowly-varying 
amplitude approximation], together with the nonlinear 
coupling coefficient from Eq. H25|l allows one to study 
nonlinear quantum optical processes such as paramet- 
ric down-conversion in the presence of realistic dielec- 
tric materials. The main advantage of our approach is 
that it automatically takes absorption — via the complex 
permittivity — and geometric boundaries — via the dyadic 
Green function — into account. The procedure to general- 
ize the theory presented above is by no means restricted 
to quadratic responses. In fact, one can construct a hi- 
erarchy of Hamiltonians with increasing number of the 
dynamical variables f (r, w) and f^(r, uj) corresponding to 
higher-order nonlinear responses. The construction en- 
sures that the equal-time commutation relations between 
the relevant field operators are preserved. We believe 
this theory represents an important step towards further 
studies with the aim to understand the ultimate limits 
on the performance of quantum optical processes. 

This work was funded by the UK Engineering and 
Physical Sciences Research Council (EPSRC). The au- 
thors thank A. Tip for helpful discussions. 
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